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Abstract 

Perturbative relations between pole and running heavy quark masses, defined 
in the Minkowski regions, are considered. Special attention is paid to the 
appearance of the kinematic vr^-effects, which exist in the coefficients of these 
series. The estimates of order QCD corrections are presented. 
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o ; 1 Introduction 

Among important parameters of QCD are the masses of c, h and t-quarks, 
^ . which are more heavy than A^l=3,4,5 number of lighter ones. They can be 

\ defined either as the poles of the corresponding renormalised heavy-quark 

^ \ propagators at = M^^^^_,_^j in the Minkowski space-like region or as the 

\ running masses m(Ar^+i)(/i^) in the MS-scheme. Their scale-dependence is 

described by the solution of the following equation 

O 



— = 



dx 



(1) 



^ \ where ^^(s) = as{s)/TT and as{s) is the QCD coupling constant of the MS- 

d \ scheme, fixed in the Minkowski reference point s > rn'^j^^_^i^, and the renor- 

malization group functions 7m(^r^+i)(a^) and f3{x) are defined as 

Pias) = '^-^ = -T.mL)ar- (3) 

The coefficients Pi^Ni) and '~ii{NL) (apart of the coefficient 70) depend from 
Nl + 1 number of active fiavours. Note, that for the MS-scheme heavy quarks 
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running masses m(iVi+i)(yU^) the Minkowskian normalization point yU^ = 
t^Inl+i) frequently used (see, e.g., Ref. [T]). In this case the definition of 
m(7v^+i) (jnf^j^ ^j^i^) may be geometrically illustrated by finding the intersection 
of the curve, which represent the inverse logarithmic scale-dependence of the 
squared running mass, with the bisectrix of the angle, formed by positive 
axises < m^^^j^^_^^-^ < oo and < /x^ < oo Q. The relations between pole and 
running heavy quark masses we will be interested read 

4 

^(^L+l) =^(WL+l)(^(7Vi+l)) Il^f (^LX(^(iVi+l)) • (4) 

n=0 

Note, that in the process of comparison of theoretical predictions for the 
e^e~-annihilation Euclidean time-like characteristic, namely Adler D-function, 
with its experimental-motivated behaviour [2j heavy quark pole masses were 
defined in the MOM-scheme, while heavy quark running masses were defined 
at the Euclidean scale fi'^ = Q^. The similar mixed MOM- MS-scheme pre- 
scriptions are also widely used to analyse heavy-quark mass dependent effects 
in characteristics of deep inelastic scattering (see e.g. [3], However, the 
processes, which may be observed at LHC, are described by theoretical pre- 
dictions in the time-like region of energies. In view of this it is important to 
study relations between different most commonly used definitions of heavy 
quark masses and to derive the relations between pole and running heavy 
quark masses, tied to the Euclidean and Minkowski regions of momentum 
transfered. This problem was analysed in Ref. [5J with the help of the special 
Kallen-Lehman type representation. Here we will consider this approach in 
more detail, presenting additional arguments in favour of theoretical back- 
ground of the investigations, performed in the work mentioned above. We 
will also update estimates of the order O(a^) terms in the relation of Eq.(jlj), 
which were obtained in Ref.|5j using the extended to the mass- dependent 
case effective-charges inspired massless approach, elaborated in Ref. 

2 Comments on application of the dispersion 
relations 

Let us discuss the subject of applicability of the Kallen-Lehman type spectral 
representations within the context of perturbative QCD. The well-defined 

^We are grateful to G.B. Pivovarov for the discussion of this topic. 
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dispersion relation for the e+e -annihilation Adler function is well known 

where UviQ"^) is the photon vacuum polarization function and R{s) ~ ImUy. 
The two-point function of the scalar quark currents rrii^N^+i^ipgipg has the 
imaginary part, which defines the scalar Higgs boson decay width into quark- 
antiquark pairs. In this case it is possible to write-down the following repre- 
sentation [7j: 

which faces no problems in the region where the the asymptotic freedom 
property of QCD holds. The same equation was used in Ref. [5J to extend the 
massless procedure of the estimates of higher-order perturbative corrections 
to the Euclidean quantities [6] to the case of Eq.([2]), which contains the 
dependence from the square of running mass m(7v^+i)((5^) defined in the 
Euclidean region. However, as was shown in Ref. [8] , the dispersive relation 
of Eq.([6]) is valid within perturbative sector only and can not be proved 
on the level of rigour, considered in Ref. [9j. Indeed, it was shown in Ref. 
[H] that in the low-energy region Eq.([n]) is ill-defined and contains fictitious 
Ag^^/Q^-term. It reflects the failure to remove the infinities from n5'(0). 
The well-defined dispersive relation, which do not contain this term, can be 
written down through the second derivative of the scalar correlator [10] . It 
leads to the following Euclidean function 

Note, however, that its perturbative expansion differs from the one, which 
corresponds to the Euclidean part of perturbative series for T{H^ — > qq), 
generated by the ill-defined in non-perturbative sector expression of Eq.([H]). 
Moreover, the application of the "approximate" dispersion relation from 
Eq.([6]) fixes the kinematic vr^-contributions to the coefficients of the per- 
turbative series for r(if° — )■ qq) both in the expanded [5] and summed up 
[8], [11] forms. Note, that the idea of the summation of vr^-terms at lowest 
order of QCD was proposed and used over thirty five years ago in the works 
of Refs. [12], [13], [Z]. 



dO^^ 
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3 Dispersion relations for the pole and 
running heavy quark masses 

Consider now the following "approximate" dispersion model of Ref . [5] for the 
heavy quark pole masses 



M, 



(Nl + I) 



i(jVi + i){m2^^^jj)+i£ 



27ri J-m(iv^+i)(m 



ds' 



■2 



T(s] 



ds 



{s + s'y 

^4 +M -n 



(8) 



with the spectral density defined as T(s) = m(jv^+i)(s) Z]n=o ""(-s) • K can 
be obtained from the dispersion-type expression for the Euclidean series 



n=0 



T{s) 



{s + Q 



2\2 



ds (9) 



where m(jv^+i)((5^) and ^^(Q^) are the heavy quark masses and the QCD cou- 
pling constant which are "running" in the Euclidean region. The applica- 
tion of Eq.(|H]) allows one to fix the relations between coefficients fni^i) and 
t^^lNi) of the perturbative series in the time-like and space-like regions as 
/<f = C Jf = tf , fiiNL) = tf{NL)+e,{NL). fiiN^) = tr(iV^)+e3(iV^), 
ffiNi) = t^^iNi) + e^lNi). The kinematic vr^-terms enter the derived in 
Ref.[S] explicit expressions for the ei(A^L)-contributions, namely 



e2(iVi 



5.89435 
3 

105.622 



-r 7oj 
0.274156A^L 



esiN^) = y |tf(/3o + 7o)(/3o + 



To 
2 



+ t 



M 



/3i7o 



+ 7i(/3o + 7o) 



10.0448A^L + 0.198001A^ 



e4(iVi 



TT 



70^ 



6 



+ 



2 

71 



y(^/3o + 7o) + y(2/3o + 7o) 



(10) 

.(11) 

(12) 
(13) 



71 



+ 72 



A) To 
2 3 



+ -^C7o(/3o + 7o)(/3o + |)(/3o + |) 

= 2272.02 - 403.951A^L + 20.6768iV^ - 0.315898Ar2 

Their A^^-dependence result from A^'^.-dependence of the coefficients Pi{Ni) 
with i > in Eq.dS}, 7i(A^L) with i > 1 in Eq.p and in Eq.dU, which 
has the following numerical form [E 

.M 



tf = 13.44396 - 1.041367Ari 



(14) 
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and comes from the analytical expression of Ref.[T2], confirmed by the in- 
dependent calculations of Ref.[Tn]. Notice, that the results of Refs.fT^. [TB] 
contain the explicit dependence from (^2 = 7r^/6-terms. The discussions 
presented above clarify that the part of these vr^-terms, explicitly visible in the 
formulae of Refs.[T5],[16], appear from the analytical continuation effect of 
Eq. fllOp . This our claim can be generalised to the level of t^-corrections, eval- 
uated analytically in Ref.[Tl] and semi-analytically in Ref.[17]. In this case 
kinematic vr^-contributions are determined by Eq. ffTT]) . The coefficients of 
the relation between heavy quark Euclidean masses, defined in the MOM on- 
shell, and MS-scheme masses, contain only remaining transcendental terms, 
typical to the on-shell scheme calculations. 

4 Estimates of corrections 

We consider now two perturbative series, namely the one of Eq. (jl]) and the 
related to it relation 

M(^,+i) = m(^,+i)(Mjv^+i)) Y: t^f (iVL)<(Mgv^+i)) . (15) 

n=0 

Keeping in mind that for < n < 3 the values of the explicit dependence 
from Nl of the coefficients iff^N^) and v^^{Ni) is already known [H], [T], we 
will study the problem of estimates of the and coefficients, using the 
effective-charges (ECH) inspired approach, developed and used in Refs.[6], 
[sj^ . It is known that the applications of this approach in the Euclidean re- 
gion at the level of and corrections give correct correct in signs and in 
order of magnitude estimates of the perturbative contributions to the num- 
ber of physical quantities (see e.g. [6], [5], [22], [23]). As to the apphcation 
of this procedure to the Minkowskian quantities, two ways are possible. The 
first, prescribes to apply the procedure of estimates in the Euclidean region 
and add explicitly calculable kinematic vr^-terms afterwards. Within the 
second way one may use the procedure of estimates in the Minkowski re- 
gion directly. It should be noted, that both ways are leading to reasonable 
predictions of signs and numerical values of perturbative series for physical 
quantities. Moreover, in the case of direct application of this approach in 
the Minkowski region, the order estimates are sometimes even closer to 
the results of the explicit calculations (see e.g. Ref. [22]). However, in the 

^ The method of ECH was proposed and developed in Refs.[TH],[in] and independently 
in Ref. [20] (see also Ref. ^l,)- 
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latter case the estimates do not reproduce the known values of the analytical 
continuation effects, similar to the ones of Eq.f lTT]) and Eq. flT^ . Note, that 
their precise knowledge is important for applying different approaches of re- 
summations of these contributions (see e.g. [2l]-|27], [8], [H]). Following two 
ways mentioned above we first estimate the values of t^lNi) coefficients and 
compare them with the results for t'^°''^^{Ni) obtained in Refs.[T3], [T7]. Sat- 
isfied by this comparison we are going one step further and estimate t^iN^)- 
coefficients, taking into account the numerical expressions for tl^°''^^{NL)- The 
concrete numbers are presented in Table 1. One can see, that the estimates 



Table 1: The estimates for tf (A^l), (jV^). 



Nl 


-/.exact 


4.ECH 


4.ECH direct 
''3 


+ECH 


J.ECH direct 


5 


73.6366 


58.0645 


48.4906 


719.339 


710.016 


4 


94.4175 


100.74 


78.243 


986.097 


1045.5 


3 


116.504 


147.303 


111.315 


1281.05 


1438.75 



obtained give correct signs and order of magnitude estimates for the values 
of (iVi)-terms. Thus, one may hope that the estimates for t^^ (Nl) are not 
far from reality. We present now concrete numbers for the coefficients of the 
series of Eq.(jl]), where for the a^-coefficients we use the estimates tf^'^^N^) 
from Table 1: 



Mr, 



mh{ml] 



1 + -as{ml) + 10.3al{ml) + 116.5a^(m2) + 1281a^(m2 
1 + ^asiml) + 9.28al(ml) + 94.4a^(m^) + 986a^(mg) 



1 + -aJm'^) + SMalimf) 



73.6a^(m2) + 719af (m^ 



(16) 
(17) 



The similar relations for Eq. ffTSj) with on-shell normalizations of running 
parameters read 



Mt^mtiMf) 



1 + -as{Ml) + l?,al{Ml) + lbQal{Ml) + 1853a^(M2) 



1 + -as{Ml) + l2al{Ml) + l?,lal{Ml) + 1460a^(M2) 
o 

1 + ^a,(M2) + llal{M^) + mal{M^) + 1101a^(M2) 



(19) 

(20) 
(21) 



The results presented in Eq.f l2CT]) give the following ratios of the squares of 
running and pole 6-quark masses 



^ 1 - -a,{M^) - 18.5559ai{M^) - 175.797a^(M,^) - 1684a^(M,^) 

(22) 

where the last term is fixed by the result of application of the ECH-motivated 
approach with adding kinematic vr^-contributions at the final step. In the 
case when the Euclidean and kinematic vr^ corrections are summed up at the 
intermediate steps, the last coefficient in Eq. fl22]) should be changed from 
-1684 to -1835. Note, that in the process of analysing the uncertainties of 
QCD predictions for r{H^ — >■ bb), preformed in the work of Ref.[28], we used 
slightly lower estimate, namely -1892. The difference is explained in part by 
smaller number of significant digits taken into account in the values of coef- 
ficients, which enter in the procedure of corresponding estimates. However, 
this difference between the values of estimated order 0{a^) contributions are 
not so numerically important. Other possible physical applications, like the 
comparison with the renormalon-based analysis of asymptotic behaviour of 
perturbative series in Eq.f lT9|) -Eq. fl2T]) (for the related theoretical discussions 
one can see Refs. [29]- [31]) are beyond the scope of this study. 
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